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ABSTRACT 

Aims. We study the effect of wind anisotropies on the stellar evolution leading to collapsars. 

Methods. Rotating models of a 60 M© star with Q/Q C rit = 0.75 on the ZAMS, accounting for shellular rotation and 
a magnetic field, with and without wind anisotropies, are computed at Z=0.002 until the end of the core He-burning 
phase. 

Results. Only the models accounting for the effects of the wind anisotropies retain enough angular momentum in their 
core to produce a Gamma Ray Burst (GRB). The chemical composition is such that a type Ic supernova event occurs. 
Conclusions. Wind anisotropies appear to be a key physical ingredient in the scenario leading to long GRBs. 
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1. Introduction 

Gamma Ray Bursts (GRB) are extraordinary events which 
can be seen at very high redshifts. The record holder is cur- 
rently the GRB050904 at z=6.29 (Cusumano et al. [2TjQ5|) . 
The study of the conditions for their appearance and the 
underlying physics may thus give important information on 
' the early Universe. 

In recent years many observations have indicated that 
long soft Gamma Ray Bursts (GRBs) are associated with 
the explosion of massive stars (see the review by Woosley 
& Bloom 2006). In a few cases it has been possible to ob- 
serve the type of the supernova associated with the GRB 
event. In four well established cases (see the review by 
Delia Valle [2006)) . the supernova was of type Ic (Galama 
et al. rreggl Mazzali et al. l2"UU6]) . This means that the spec- 
tra do not show any hydrogen and helium lines. Such ob- 
servations gave strong support to the collapsar model pro- 
posed by Woosley (|1993p . In this model, the progenitor of 
the GRB is a massive Wolf-Rayet (WR) star that collapses 
into a fast rotating black hole. Fast rotation is important 
since it allows part of the infalling matter to form a disk 
around the black hole. The disk acts as an efficient engine to 
extract gravitational energy. Part of this energy is used to 
power strong polar jets. The GRB arises from interactions 
in these jets. The short duration of the GRBs points toward 
a compact progenitor, i.e. a star having lost its outer enve- 
lope, a WR star possibly of the WO subtype (Woosley 1993; 
Mazzali et al. |2"UTJ61) . 

While this general scenario is now well accepted, the 
precise features of the long GRB progenitors are still un- 
known. First, not all WR stars will give birth to a GRB. 
Delia Valle (J2006) for instance estimates that the number 
ratio of GRBs to type Ibc SNe is between 0.5-4%. Hammer 
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et al. (|2006|) suggest that GRB progenitors are runaway 
massive stars ejected from compact massive star clusters. 
Fruchter et al. (2006) find that GRBs do not occur in similar 
environments to supernovae: they found that the GRBs are 
far more concentrated in the very brightest regions of their 
hosts galaxies than are supernovae. Moreover, the GRB 
host galaxies are fainter and more irregular than the hosts 
of supernovae. These observations may point toward a very 
massive progenitor in metal poor environments. In seven 
cases the metallicity of the galaxy hosting the GRB was 
measured or at least was given an upper limit. In all the 
cases, the metallicity was found inferior to one third solar 
(reported as a personal communication from Vreeswijk in 
Fruchter et al. [20061) . 

The progenitors of long GRBs represent a real challenge 
for stellar physics. On one hand, it is necessary that the star 
loses a large fraction of its initial mass in order to make a 
WO star. On the other hand, the star must keep a high spe- 
cific angular momentum j (> 10 16 cm~ 2 s -1 ) in its core for 
the collapsar model to work. These two constraints are rela- 
tively contradictory, because removing a lot of mass implies 
losing a lot of angular momentum. This is the difficulty for 
models of GRB progenitors. 

Models accounting for the effects of rotation but without 
magnetic fields have no difficulty in preserving enough an- 
gular momentum in their central regions (Heger et al. 2000; 
Hirschi et al. I2005j) . The conditions of having WO stars as 
progenitors is consistent with metal poor environments (see 
Hirschi et al. 2005). However, these models may somewhat 
overestimate the frequency of GRBs and cannot reproduce 
the observed rotation rate of young pulsars. When the mag- 
netic fields are accounted for (Heger et al. 2005) as proposed 
in the dynamo theory of Spruit (|2002p . the formation of a 
collapsar is difficult, at least for not too extreme velocities. 
In that case, the magnetic field imposes solid body rotation 
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Fig. 1. Evolutionary tracks for rotating models with and 
without wind anisotropies. The plotted effective tempera- 
tures correspond to the surface averaged effective tempera- 
ture at the hydrostatic surface (no correction for the opti- 
cal thickness of the wind) . The beginning of various phases 
are indicated along the evolutionary sequences: Bup, WR, 
and X c = indicate when the star encounters the critical 
or break-up limit, enters the WR phase and when core H- 
burning stops. The core He-burning phase begins at the 
point labeled He-burning. The transition phase between 
the WNL and WC phases called the WN/WC phase corre- 
sponds to the portion of the tracks comprised between the 
two tick marks with, in the middle, the label WN/WC. 
After that point the star is a WC star. The entrance into 
the WO phase is indicated (if it occurs), as well as the 
endpoint of the core He-burning phase labeled by Y c = 0. 



during the Main Sequence phase, allowing efficient transfer 
and extraction of angular momentum from the core and the 
models have not enough angular momentum in their core 
to give birth to a collapsar. 

A way to overcome this difficulty was proposed by Yoon 
& Langer (j2"Ut?5|) and Woosley & Heger (|2T)U6T) . who com- 
puted the evolution of very fast rotating massive stars with 
homogeneous evolution similar to the models computed by 
Maeder (1987). These stars remain in the blue part of the 
HR diagram and avoid the red supergiant phase. However, 
such models enter quite early into the WR phase due to the 
action of strong mixing and thus have to escape the danger 
of losing too much angular momentum in the very strong 
WR stellar winds. To avoid this, a second crucial clement 
has to be accounted for: the fact recently shown that the 
winds of WR stars are metallicity dependent (Crowther et 
al. 120021 Vink & Koter |2"UT)5)) . i.e. weaker at lower metal- 
licity. When such mass loss recipes are used, the homoge- 
neous models may lead to collapsars. Within this scenario, 
the rarity of the GRB events is due to the fact that the 
GRB progenitors only occur for high initial stellar masses, 
for very high rotational velocities and for low metallicities. 



Models in the literature do not account for the effects of 
the wind anisotropies induced by rotation. Here we present 
the first computation including magnetic field and wind 
anisotropies. We revisit the above scenario (homogeneous 
evolution) accounting for this effect. In Sect. 2, the model 
ingredients are presented. In Sect. 3 to 5 we discuss the 
effects of wind anisotropies on the evolution, rotation and 
chemical abundances. Conclusions are given in Sect. 6. 

2. The physical ingredients 

The stellar ingredients are the same as in Meynet & 
Maeder (2005|) with the following exceptions: 1) For the 
mass loss rates we use the same prescriptions as in Yoon & 
Langer (2005; their WR1 recipe during the WR phase). We 
suppose that the star enters the WR phase when the hy- 
drogen surface abundance becomes inferior to 0.3 in mass 
fraction; 2) The effects of the magnetic field are accounted 
for as in Maeder & Meynet (|2005D . Let us recall that in 
rotating models with magnetic field, the main mechanism 
driving the transport of the chemical species is meridional 
circulation. 

We calculate two 60 M Q models with an initial composi- 
tion X = 0.748 and Z = 0.002 with a value n/n clit = 0.75 
on the ZAMS (v ini = 523 km s _1 ). Note that such high 
initial velocities are required for at least two purposes: first 
in order to have a well developed wind anisotropy, second 
such initial velocities seem to be required to obtain a col- 
lapsar in rotating models with magnetic fields (see Yoon 
& Langer I21K351 Woosley & Heger |2"0"01)]) . Faster rotation 
will likely give similar results to those shown here, while 
lower initial rotation would probably show no important 
wind anisotropy effect and thus would not lead to a col- 
lapsar. One model was computed assuming isotropic winds 
and one model with anisotropic winds. The effects of the 
wind anisotropies are computed as in Maeder (120021) . The 
computations were followed until the end of the core He- 
burning phase. After that stage, the angular momentum 
of the core does not change much (see Hirschi et al. 2005) 
and has already reached a value near the one it will have 
at the presupernova stage. This conclusion was based on 
models without a magnetic field and with a lower initial 
rotation. Woosley & Heger (|2006|) . however show that the 
specific angular momentum at the Lagrangian mass of 3M Q 
decreases by an order of magnitude between the end of the 
core He-burning phase and the presupernova stage in mod- 
els with a magnetic field and fast rotation. In similar models 
by Yoon and Langer (|2005p . on the contrary, the decrease 
of the inner specific angular momentum between the end of 
the He-burning and the end of the core C-burning is very 
modest (cf. their Fig. 5). Thus, both our models (Hirschi et 
al. l2005p and those of Yoon & Langer (2005 ) support a small 
change of the specific angular momentum in the advanced 
nuclear stages, in contrast with Woosley & Heger (2006). 
Theoretically, we do not see which mechanism (among those 
included) could efficiently remove the specific angular mo- 
mentum from the core in the late nuclear stages, since the 
two most efficient processes (i.e. magnetic field and trans- 
port by meridian circulation) vanish or become inefficient 
in regions with ^-gradients as strong as those near the core 
of massive stars. 

As will be discussed below, the surface velocity of our 
"wind anisotropic" model reaches the critical limit. At the 
critical limit, wind anisotropies are not accounted for and 
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Fig. 2. Evolution as a function of time of the total masses 
(upper curve) and the masses of the convective cores (lower 
curve) for the models with isotropic winds (upper panel) 
and anisotropic winds (lower panel). The isotropic model 
reaches X c = at 4.94 x 10 6 yr and Y c = at 5.29 x 10 6 yr. 
For the anisotropic case, the corresponding ages are 5.29 x 
10 6 yr and 5.78 x 10 6 yr. The labels Bup, WR and End 
Bup show the points when the star encounters the Break- 
up limit, enters the WR phase and evolves away from the 
Break-up limit. 



the loss of angular momentum corresponds to that the star 
would have from a spherical wind. For a star of T e ff < 25000 
K, this would probably be wrong, because the opacity at 
the equator is higher due to the lower T e g and it would 
drive the matter outward in the equatorial plane. For very 
hot stars with T e fi > 50000 K, the star keeps a strong polar 
wind at the critical limit, and some matter is ejected in the 
equatorial plane. It is difficult at the present time to predict 
what would be the net effect on the angular momentum 
content of the star. However, at the critical limit, it is in 
general sufficient that a small amount of mass is removed 
to bring the star away from the critical limit, i.e. within the 
regions of polar stellar winds. Therefore the combination of 
anisotropic winds and spherical winds at break-up may be 
an acceptable assumption, at least as long as the physics of 
break-up has not been more investigated. 



Fig. 3. Evolution as a function of time of the equatorial ve- 
locity (top) and of the ratio 0/r2 cr ;t at the surface (middle). 
The bottom panel shows the evolution of the specific angu- 
lar momentum at the Lagrangian mass coordinate M r = 3 
M . 

3. Effects of wind anisotropies on the evolution 

Since homogeneous evolution may have observational coun- 
terparts (see e.g. Walborn et al. I2004p . we discuss some 
caracteristics of our evolutionary tracks in the HR diagram 
(see Fig. [1]). Due to their fast rotation and thus strong in- 
ternal mixing, the stars remain in the blue part of the HR 
diagram. Both models follow qualitatively the same evolu- 
tion, going from an O-type star phase through a WR phase, 
although with some striking differences: 

1) The "anisotropic" model encounters the break-up or 
critical limit during the MS phase. This is of course a nat- 
ural consequence of the smaller amount of angular momen- 
tum lost. 

2) The entrance into the WR phase occurs at about the 
same evolutionary stage in both models, i.e. with a mass 
fraction of hydrogen at the center X c — 0.28 and X c = 0.29 
in the "isotropic" respectively "anisotropic" wind stellar 
models and with an actual mass of 52.4 M Q , respectively 
55.7 Mq. However in the case of the anisotropic model, the 
first part of the WR phase (between the labels WR and 
X c = in Fig. []} shows a steep decrease in luminosity, 
while in the isotropic model the corresponding change is 
limited. The reason is that in the anisotropic model dur- 
ing this phase, the star continues to be at the critical limit 
(see Fig. [3]), where it loses more mass than the isotropic 
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Fig. 4. Abundances (in mass fractions) as a function of the 
Lagrangian mass of various elements inside stellar models 
without (top) and with (bottom) wind anisotropies when 
F r = 0. 



model which is far from this limit. Thus, due to the mass- 
luminosity relation followed by homogeneous models, the 
"anisotropic" model undergoes a more pronounced lumi- 
nosity decrease. When X c = 0, the isotropic model has an 
actual mass of 35.1 M Q , while the anisotropic model has 
only 22.1 M Q (see also Fig. EJ. 

3) The transition between the end of the core H-burning 
and the beginning of the core He-burning phase is different 
in the two models (cf. Fig.[T]). This comes from the fact that 
the isotropic model still has a small H-rich envelope (mass 
fraction of hydrogen at the surface X s = 0.07), while the 
anisotropic model is nearly a pure He-core (X s = 0.001). 
When, at the end of the core H-burning phase, a small H- 
rich envelope is present, a phase of expansion follows the 
contraction of the star. The expansion phase is stopped 
when the star encounters the f2r-limit. This corresponds 
to the point labeled "He-burning" along the isotropic track 
in Fig. [T] This produces an enhancement of the mass loss 
for a short time. After this phase, the star describes some 
back and forth movements in the HR diagram, which we 
have schematically represented by a hatched zone in Fig. [1] 
For the anisotropic model, the core simply contracts. More 
mass is lost since contraction bring more mass above the 
critical limit. This increase in the mass loss finally makes 
the star evolve definitively away from the critical limit. 



4) In both models, the end of the core He-burning phase 
occurs away from the critical limit (Fig. [3]). The evolutions 
of the two models (after the label WN/WC) differ mainly 
in the luminosity: the anisotropic model has a luminosity 
about 0.30 dex below that of the isotropic model, because 
the first model has lost much more mass than the second 
one. The final masses of the isotropic and anisotropic mod- 
els are 19.3 and 11.5 M© respectively. 



4. Evolution of rotation 

Figure [3] shows the evolution as a function of time of the 
equatorial velocity at the surface, of the ratio of the sur- 
face angular velocity to the critical velocity and of the spe- 
cific angular momentum at the Lagrangian mass coordinate 
M r = 3 Mq. In the anisotropic model, the surface velocity 
increases first and then remains at the critical limit during 
the rest of the core H-burning phase. When the star enters 
the WR phase, the star continues to encounter the critical 
limit, however the higher mass loss rates maintain the star 
on average at a slightly greater distance from the critical 
limit. In the isotropic model, the surface velocity decreases 
and approaches the critical limit only during a very short 
episode in the transition phase between the core H- and 
He-burning phase. In the bottom panel of Fig. [3l one sees 
the importance of the wind anisotropies on the specific an- 
gular momentum at M r = 3 Mq. As long as the stellar 
models do not lose too much mass by stellar winds, i.e. 
before the WR phase, the specific angular momentum at 
M r = 3 Mq remains nearly constant. When the star enters 
the WR phase, the specific angular momentum strongly 
decreases in the isotropic model, while, when anisotropic 
winds are considered, the specific angular momentum shows 
nearly no decrease. This allows the star to enter the core 
He-burning phase with a high angular momentum content. 
At M r = 3 Mq, the specific angular momentum is about 9 
times larger in the anisotropic model than in the isotropic 
one. Moreover, since the anisotropic model has lost more 
mass (but much less angular momentum), the mass loss 
rates in the core He-burning phase (due to their specific de- 
pendence on mass) are lower than for the isotropic model. 
This allows the anisotropic model to end its evolution with 
sufficient angular momentum to give birth to a collapsar. 
On the contrary the isotropic model has too low a central 
specific angular momentum to be a progenitor of a long 
GRB. 



5. Chemical abundances 

In Fig. [H the variations of the abundances in the two 60 
Mq models are shown when Y c = 0. The isotropic model 
has less He than CO at the surface, while for the anisotropic 
model it is the opposite. This at first sight might be surpris- 
ing. Indeed, why is the model losing the greatest quantity of 
mass, ending with more helium at the surface? The reason is 
that when the star loses more mass, it enters the WC phase 
at an earlier stage of the core He-burning phase, i.e. when 
smaller amounts of helium have been processed into carbon 
and oxygen (Smith & Maeder 4991). This is what happens 
here. The anisotropic model enters the WC phase when 
the mass fraction of helium at the center Y c = 0.93, while 
the isotropic model enters the WC phase with Y c = 0.57. 
Taken at face value, the fact that the anisotropic model 
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produce a WC progenitor differs with the type Ic super- 
nova associated with long GRB. However if we integrate 
the total quantity of helium ejected (supposing that all the 
He-rich layers are ejected at the time of the GRB event), 
we obtain not too different values, about 0.53 M Q for the 
isotropic model and 0.77 M for the anisotropic one. Such 
quantities of helium, compared to current supernovae, may 
escape detection during a type Ic event and agree with the 
idea that type Ic supernovae either eject He-deficient ma- 
terial or fail to non thermally excite their helium (Woosley 
& Eastman [TM7I)- 



Woosley, S.E., Eastman, R.G. 1997, in "Thermonuclear Supernovae", 
P. Ruiz-LaPuente, R. Canal & J. Isern eds., Dordrecht, Kluwer, 
p. 821 

Woosley, S. E., Heger, A. 2006, AIP Conference Proceedings, 836, 
p.398. 

Woosley, S.E., Bloom, J.S. 2006, ARAA, 44, 507 



6. Conclusions 

We have shown that wind anisotropies play an important 
role in the scenario leading to collapsars and then to long 
GRBs. Wind anisotropies help maintain a high content of 
angular momentum in the stellar core, while they make 
the star lose very large amounts of mass. The anisotropic 
model loses more mass than the isotropic model, because 
the anisotropic model remains near the critical limit for 
most of the core H-burning phase. Interestingly, despite 
these important losses of matter, the anisotropic model 
ends its stellar lifetime with high surface velocities, be- 
tween 200-400 km s" 1 , while the isotropic model shows 
very modest values (a few tenths of km s _1 , see Fig. 
At very low metallicity, GRB progenitors are expected to 
show high surface ratios of [Ne/Fe]. Usually the abundance 
of 22 Ne observed at the surface of a WC or WO star is 
equal to the initial CNO abundance. In the case of primary 
nitrogen production (expected to be abundantly produced 
in very metal poor fast rotating stars), the abundance of 
22 Ne might be much higher that the initial CNO abun- 
dance. Thus Wolf-Rayet progenitors of collapsars at very 
low metallicity should show, together with a high surface 
rotational velocity, a very high [Ne/Fe] ratio. In a further 
work, we will explore the consequences of wind anisotropies 
on a larger ranges of masses, metallicities and initial rota- 
tion velocities. 
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